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ABSTRACT 


This  report  describes  accomplishments  in  the  period  1977-1978  in  basic  research 
dealing  with  phenomena  in  parallel  viscous  flows.  An  experiment  has  been  com- 
pleted on  the  development  of  the  mean  flow  in  wakes  shed  by  nozzle  arrays 
typical  of  gas-dynamic  lasers.  In  this  experiment,  the  wake  velocity,  temper- 
ature and  density  was  mapped  in  wakes  generated  at  Mach  4 in  the  Company's 
Supersonic  Wind  Tunnel.  The  results  are  in  general  agreement  with  theoretical 
predictions  developed  under  an  earlier  OSR  study  in  which  theoretical  rules  on 
wake  flows  were  set  up  for  laser  cavity  design.  A second  problem  addressed 
dealt  with  a novel  approach  to  predicting  transition  in  parallel  shear  flows. 

The  basic  physics  in  this  approach  consists  of  the  postulate  of  a universal 
turbulence  Reynolds  number,  below  which  self-preserving  turbulence  is  impossible. 
Although  this  condition  is  only  one  of  necessity,  formulas  derived  with  it  had 
earlier  predicted  wake  transition  remarkably  well.  In  the  present  period,  the 
postulate  was  applied  to  flat  plate  boundary  layers  in  the  range  0 < M < 10 
with  and  without  heat  transfer.  The  computations  uncovered  a mechanism  for 
the  long  unexplained  "transition  reversal"  phenomenon  and  otherwise  conformed 
to  earlier  transition  observations.  Finally,  the  data  analysis  was  completed 
on  an  earlier  hypersonic  boundary  layer  stability  experiment  done  at  AEDC. 

The  data  show  little  change  of  the  stability  diagram  with  heat  transfer  although 
the  amplification  rates  were  substantially  increased  when  the  wall  temperature 
was  lowered.  The  present  results  were  disseminated  through  two  technical  reports, 
three  journal  articles,  three  delivered  papers  and  by  personal  contact  with  Air 
Force  engineering  centers. 
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INTRODUCTION  AND  BACKGROUND 


The  work  in  progress  in  this  laboratory  under  AFOSR  Contract  F44620-75-C-0016 
addresses  basic  problems  in  parallel  shear  flows  of  compressible  fluids.  These 
problems  are  well  defined  and  consist  of  the  issue  of  hydrodynamic  stability  of 
high-speed  boundary  layers,  the  issue  of  transition  from  the  laminar  to  the 
turbulent  state  of  a shear  flow,  and  the  issue  of  the  prediction  of  free  shear 
flow  characteristics  as  encountered  in  engineering  devices. 

In  contrast  to  many  other  OSR-supported  studies,  the  present  work  is  done  in 
an  industrial  organization  involved  in  the  engineering  of  many  modern  systems 
for  defense  and  civilian  needs.  This  enhances  the  relevance  of  the  research 
performed,  makes  the  applicability  of  its  results  highly  visible,  and  fre- 
quently allows  the  testing  of  the  research  conclusions  by  comparison  with  the 
performance  of  actual  systems.  Furthermore,  the  Fluid  Mechanics  Section  of 
this  Company  is  engaged  in  contract  research  for  military  directorates  other 
than  AFOSR;  these  projects  are  closely  related  to  the  subject  of  the  present 
contract,  and  provide  motivation,  direction  and  testing  grounds  for  the  present 
work.  Specific  examples  of  this  interaction  are  mentioned  in  this  report. 

The  present  report  covers  progress  in  the  period  from  October  1977  through 
September  1978.  The  report  addresses  separately  the  subjects  covered,  namely 
(a)  development  of  wakes  in  gasdynamic  lasers,  (b)  hydrodynamic  stability  of 
hypersonic  boundary  layers  and  (c)  a novel  theoretical  approach  to  the  laminar- 
turbulent  transition  program.  The  status  of  documentation  and  publication  of 
results  is  reported,  and  interaction  with  other  individuals  and  agencies  is 
indicated. 

2.  THE  WAKES  OF  GASDYNAMIC  LASER  NOZZLE  CUSPS 
2.1  BACKGROUND  AND  PREVIOUS  WORK 

Our  involvement  in  the  nozzle-cusp  wake  problem  was  motivated  by  the  need  to 
organize  the  knowledge  of  the  two-dimensional  wakes  trailing  the  multiple  cusps 
in  gasdynamic  laser  cavities.  For  such  lasers,  the  wakes  divide  homogeneous 
stream,  and  interest  in  die  wake  arises  mainly  because  of  possible  attenuation 
and  defocusing  of  the  beam.  Knowledge  of  compressible  wakes  was  already  adequate 
for  no  heat  addition  in  the  volume,  if  the  wake  could  be  characterized  in  the 
momentum  sense  by  a drag  integral  and  in  the  energy  sense  by  a heat  integral. 

What  was  missing,  however,  was  a consideration  of  the  entire  cusp-wake  system, 
an  accounting  of  the  sidewall  effects,  and  the  question  of  far  wake  interaction. 
The  cusp-wake  system  is  important  to  the  designer  who  needs  to  compute  wake 
development  in  terms  of  plenum,  nozzle  geometry  and  heat  transfer  conditions. 

In  the  period  1975-77,  we  began  systematically  to  review,  calculate  and  compile 
into  compact  form,  the  characteristics  of  two-dimensional  laminar  and  turbulent 
compressible  wakes  as  generated  by  flows  of  arbitrary  Mach  and  Reynolds  numbers, 
wall  temperature,  specific  heat  ratio  and  other  molecular  properties.  These 
results,  which  appeared  in  report  form  (References  1 and  2)  and  in  journal 
publication  form  (References  3 and  4),  revealed  some  unexpected  phenomena  such 


as  Che  peculiar  behavior  of  a wake  which  is  turbulent  starting  at  the  cusp 
trailing  edge  ("transitional"  wake)  and  the  irregular  variation  of  its  fluctu- 
ation magnitudes  with  Reynolds  number.  To  arrive  at  these  results,  a combination 
was  used  of  known  flow  properties  (e.g.,  laminar  and  turbulent  diffusion  rates) 
and  of  some  new,  key  assumptions;  for  example,  it  was  assumed  that  the  virtual 
wake  origin  lies  at  the  nozzle  throat  for  the  fully  laminar  and  fully  turbulent 
wake. 

To  test  the  predictions  of  References  1 through  4,  a triple  nozzle  arrangement 
was  set  up  in  the  Aeronutronic  Supersonic  Wind  Tunnel  to  produce  exit  M ■ 4 and 
nozzle  wall  temperatures  ^ down  to  0.7  T0.  The  purpose  was  to  measure  the  cusp 
wake  properties  in  this  set-up  and  compare  them  with  the  properties  (velocity, 
temperature,  etc.),  of  the  above  references.  By  late  1977,  the  wind  tunnel 
set-up  was  working  smoothly;  preliminary  results  on  the  mean  flow  had  been 
obtained  and  reported  in  Reference  5.  The  final  experiments,  which  had  been 
scheduled  for  completion  in  the  1977-78  period,  are  mentioned  below. 

2.2  NOZZLE  CUSP  WAKE  MEASUREMENTS  IN  THE  PRESENT  WORKING  PERIOD 

In  the  present  period,  the  wake  measurements  were  completed  and  are  reported  in 
Reference  6.  The  results  are  aptly  sumnarized  in  that  reference: 

"Experimental  measurements  are  described  in  the  supersonic  wakes  of  cusps 
formed  synme trie ally  by  DeLaval  nozzle  contours.  The  nozzle  exit  air  flow, 
serving  as  the  wake  edge  condition,  was  continuous  at  Mach  4 and  nominal 
Reynolds  number  of  100,000.  The  measurements  spanned  Che  three  cases  of  the 
adiabatic  laminar  wake,  the  adiabatic  laminar  wake  undergoing  turbulent  tran- 
sition downstream  of  the  trailing  edge,  and  the  laminar  wake  produced  by  a 
nozzle  cusp  cooled  to  0.7  TQ.  The  mean  flow  properties  of  such  wakes  were 
predicted  by  the  author  in  an  earlier  theoretical  study  and  the  present  work 
aimed  at  providing  a test  of  Che  theory.  Verification  was  indeed  obtained  for 
key  assumptions  in  the  theory  such  as  che  location  of  Che  virtual  origin  and 
Che  existence  of  laminar  similarity  and  turbulent  self-preservation.  As  pre- 
dicted, che  nozzle  cusp  temperature  affects  che  temperature  and  density,  but 
not  the  fluid  velocity.  The  non-equilibrium  region  associated  with  the  trailing 
edge  was  f >und  to  extend  to  about  14  boundary  layer  thicknesses.  Chief  source 
of  anomalies  was  Che  network  of  crailing-edge  source  and  their  reflections  which 
introduce  pressure  gradients  suspected  to  be  responsible  for  an  overall  de- 
pression of  the  velocity  defect.  The  disparities  are,  nevertheless,  small  and 
the  agreement  between  theory  and  experiment  is  considered  adequate." 

2.3  DISCUSSION  AND  UTILITY  OF  THE  RESULTS 

The  results  obtained  have  largely  validated  and  reinforced  the  theory  of  Ref- 
erences 1 through  4;  by  this,  the  present  author  implies  chat  he  would  un- 
hesitatingly use  Che  references  to  predict  trends  and  magnitudes  of  wake  velocity, 
temperature,  density  and  turbulence  magnitudes  in  most  gasdynamic  laser  cavities 
configured  like  che  test  of  Reference  6.  While  much  of  Che  predictions  of 
References  1,  etc.,  were  never  in  doubt  (e.g.,  Che  turbulent  diffusion  rates), 
some  findings  of  che  cescs  of  Reference  6,  such  as  that  concerning  che  virtual 
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origin,  are  new  contributions  of  practical  importance. 

Following  publication  of  Reference  6,  Peterson  (Reference  7)  completed  and 
published  results  of  a similar  experiment  at  M 6 and  at  two  different  Reynolds 
numbers.  He  found  that  the  turbulent  theory  of  Reference  1 agreed  with  his  tur- 
bulent results,  but  that  his  laminar  results  lay  somewhere  between  our  laminar 
and  transitional  predictions.  Disagreement  by  about  50%  in  the  velocity  and 
temperature  defects  between  theory  (Reference  1)  and  the  data  were  also  noticed 
by  us  in  the  experiments  of  Reference  6.  Two  possible  explanations  exist.  One 
is  that  shock-wave  interactions  in  the  multiple  nozzles  are  important;  in  fact, 
in  Reference  6,  it  is  noted  chat  the  longitudinal  pressure  along  the  wake  ex- 
periences considerable  variations.  Peterson  (Reference  7)  noticed  a similar 
event.  The  theory,  on  the  other  hand,  is  formulated  for  constant-pressure  wakes, 
while  pressure  gradients  can  be  generated  by  the  criss-crossing  shocks  and  wakes. 
The  second  possibility  is  that  the  laminar  wakes  studied  (in  both  ours  and 
Peterson's  cases)  are  in  fact,  transitional.  Support  for  this  possibility  comes 
from  the  finding  that  the  transition  theory  we  advanced  (see  Section  3)  predicts 
the  laminar  wake  investigated  in  Reference  6 to  be  on  the  threshold  of  becoming 
turbulent.  Transitional  wakes  indeed  have  lower  velocity  defects,  exactly  as 
found . 

With  the  satisfactory  outcome  of  Che  experiments  of  Reference  6,  this  OSR  effort 
has  produced  a compendium  of  wake  characteristics,  supported  by  experiments, 
which  should  become  a convenient  reference  for  the  discussion  of  the  mean- flow 
properties  of  laser  cusp  wakes.  Furthermore,  the  mean- flow  results  provided  us 
with  the  basis  of  making  estimates  of  the  turbulent  magnitudes,  which  are 
presented  in  References  2 and  4.  The  outstanding  result  of  these  turbulence 
predictions  is  that  the  rms  wideband  density  fluctuation  on  the  wake  plane  of 
symmetry  ("axis")  reaches  a maximum  of  177.  of  the  mean  density  in  Che  cavity 
regardless  of  Mach  number,  Reynolds  number,  heat  transfer,  etc.  Shortage  of 
time  prevented  extensive  turbulence  measurements  in  the  present  contract;  it 
may  even  be  argued  that  the  existing  background  in  the  turbulence  properties 
of  adiabatic  wakes  (see  Reference  1 for  example)  does  not  merit  further  work 
on  this  subject.  However,  there  is  justification  for  measuring  the  turbulence 
properties  for  cooled  nozzle  cusps,  since  no  turbulence  data  for  compressible 
wakes  with  heat  transfer  are  available.  It  is  planned  to  make  at  least  some 
such  measurements  in  the  SWT,  during  the  next  contract  period  using  the  hard- 
ware described  in  Reference  6. 

The  double  wake  system  with  which  the  experiment  of  Reference  6 was  carried 
out  in  Aeronutronic's  Supersonic  Wind  Tunnel  is  pictured  in  Figure  1.  Typical 
results  of  this  experiment  are  shown  in  Figure  2 in  which  the  theoretical  pre- 
dictions of  Reference  6 are  also  drawn. 

3.  THE  DISSIPATION  CRITERION  OF  TRANSITION 

3.1  INTRODUCTION 

After  twenty  years'  involvement  in  transition  research,  the  present  author  has 
become  skeptical,  if  not  pessimistic,  about  the  prospects  of  a "theory"  of 
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transition,  mainly  due  to  the  following  observations. 

First,  there  is  a persistent  and  large  gap  between  the  physical  and  mathematical 
approaches  to  the  transition  problem  on  one  hand,  and  the  practice  of  transition 
prediction  on  the  other.  Fluid  dynamicists  have  been  postulating  theoretical 
avenues  and  performing  carefully  controlled  experiments,  in  seeming  disregard 
of  the  engineers  and  designers  who  invent  their  own  rules  for  predicting  tran- 
sition in  practical  situations  - and  vice  versa. 

Secondly,  the  pace  of  the  theoretical  and  microscopic  experimental  work  on  the 
transition  problem  is  inadequate  to  handle  the  large  number  of  critical  design 
problems  for  which  transition  prediction  is  an  everyday  necessity.  For  ex- 
ample, boundary  layer  stability  analysis  and  experiment*  have  been  active  for 
over  forty  years  and  yet  no  transition  prediction  from  the  stability  findings 
alone  is  yet  possible.  It  is  indeed  remarkable  that  aircraft  and  missiles 
have  been  successfully  operating  in  a spectrum  of  flight  parameters  (Mach  and 
Reynolds  numbers,  heat  transfer,  etc.),  for  so  long  without  the  benefit  of  a 
rational  solution  to  the  transition  problem. 

Third,  many  current  predictions  of  laminar- turbulent  transition  quoted  by 
engineers,  are  invariably  based  on  empirical  formulas  which  are  ephemeral  in 
their  history  and  severely  limited  in  their  scope.  Typical  are  cases  where 
these  formulas  derive  from  curve- fitting  of  data  from  indirect  measurements  of 
dubious  validity  taken  in  uncontrolled  experiments.  Such  approaches  are  rarely 
supported  by  dimensional  analysis  and  feature  mysterious  numerical  constants 
which  require  adjustment  every  time  the  flow  parameters  or  regimes  change. 

These  thoughts  prompted  the  author  to  search  for  a method  of  predicting  the 
onset  of  transition  to  turbulence  which  is: 

(a)  based  on  accepted  physics, 

(b)  mathematically  simple, 

(c)  valid  for  all  coiranon  forms  of  shear  flows, 

(d)  compatible  with  past  knowledge  of  shear  flow  behavior,  i.e.,  utilizing 
established  observations  of  the  transition  behavior  and  the  character  of 
laminar  and  turbulent  flows, 

(e)  independent  of  the  microscopic  mechanism  of  transition. 

With  these  provisions  one  is  led  almost  naturally  to  seek  a minimal  condition, 
or  threshold,  for  the  existence  of  turbulence  in  a flow.  We  are  thus  forced  to 
think  in  elementary  terms  and  one  such  elementary  path  is  to  consider  the  mini- 
mum possible  width  of  a turbulent  zone.  Since  such  a zone  contains  secondary 
motions  of  a continuum  fluid,  it  is  clearly  unable  to  reach  the  small  dimensions 
possible  of  a laminar  fluid;  the  factor  controlling  the  width  of  such  a zone, 

*Stability  experiments,  one  of  which  is  done  under  the  present  program  (Section 
4),  are  invaluable  as  the  only  rational  approach  to  finding  the  mechanism  of 
transition.  What  is  stressed  here  is  the  need  for  complementary  approach  which 
can  predict  the  occurrence  of  transition  with  no  sacrifice  of  physical  insight 
and  rigor. 
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or  Che  size  of  eddies,  must  be  connected  with  the  molecular  or  viscous  effects. 
One  is  thus  led  to  the  familiar  concept  of  the  viscous  dissipation  threshold, 
below  which  turbulent  motion  is  forbidden  by  such  dissipation.  In  abbreviated 
thinking,  this  recalls  the  equally  familiar  turbulence  Reynolds  number 


ReA 


(1) 


where  u'  is  the  velocity  fluctuation,  A the  correlation  length  scale  and  v the 
kinematic  viscosity.  The  limit  we  seek  is  a limit  to  this  number. 

3.2  VARIATION  OF  THE  TURBULENCE  REYNOLDS  NUMBER 

The  quantity  Re/y  will  be  the  focus  of  our  discussion  and  it  will  be  useful  if 
we  first  consider  its  magnitude  and  development  generally,  as  affected  by  (and 
as  affecting)  the  common  turbulent  shear  flows.  By  definition,  such  flows 
feature  some  type  of  balance  between  gain  and  loss  of  turbulence  and  this 
balance  allows  the  scaling  of  flow  properties  by  similarity  laws.  Furthermore, 
we  will  temporarily  consider  isobaric.  flows  (J3p  3 0)  since  the  extension  to 
pressure  gradients  makes  the  initial  discussion  cumbersome  without  adding 
insight.  The  scaling  and  the  isobaric  assumption  simplify  the  first  overview 
of  Che  behavior  of  Re^;  later  embelishments  will  not  change  the  picture. 

We  can  gain  considerable  intuition  by  first  developing  an  overview  of  the 
variation  of  Re.  in  the  common  turbulent  shear  flows  such  as  the  wake,  the 
mixing  layer,  the  jet  and  the  boundary  layer.  Much  of  this  intuition  can  be 
obtained  by  first  neglecting  the  details  of  transverse  variations  of  u‘, 

A and  v.  For  the  moment  let  us  therefore  utilize  only  the  basic  features  of 
turbulent  similarity,  that  u1  and  A bear  certain  relations  to  the  velocity 
scale  u(0)  - ue  and  the  width  scale  b;  here  u(0),  ue  and  b are  the  lowest  and 
highest  mean  velocity  in  the  flow  and  the  flow  width,  respectively.  We  can 
then  write 


ue-u(0) 


u ' 3 Ct  [ue  - u(0)]  = Ct  ( — - ) ue. 


C.  = const 


(2) 


A = r b ■ C„b  , C-  * const 


(3) 


In  the  same  spirit,  we  can  write 

n 


v - — ve  - ( ~)  v , n > 0 

v o ^ e 


(4) 


Le 
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The  net  result  is 


u -u(0)  T-n  T "n 

ReA  ” C3  ( — Z ^ b ^Reb^  ^ f”  ) * C3  wb  Re'  ( ^ ),  » const 

(5) 

where  w s (ue  - u(0))/ug  is  the  velocity  defect  of  the  flow,  and  Re'  the  unit 
Reynolds  number  which  is  a fixed  number  for  each  given  flow. 

This  simple  formula  is  already  sufficient  to  give  us  a general  picture  of  the 
turbulence  Reynolds  number  variation  in  our  typical  flows.  The  product  wb  is 
dependent  on  x,  i.e.,  the  distance  from  the  origin  of  the  flow  (e.g.,  the 
exit  of  a jet,  or  the  wake-shedding  body);  it  depends  on  the  axisymmetric  or 
two-dimensional  geometry  of  the  flow.  The  Re*  fixes  the  characteristic 
Reynolds  number  of  the  problem,  and  the  temperature  ratio  contains  the  com- 
pressibility and  heat  transfer  information.  Thus  Re^  is  found  to  depend  on 
distance,  on  geometry,  on  Reynolds  number  and  on  compressibility  and  neat 
transfer. 

In  incompressible  flows,  defined  as  T * Te  in  Equation  (5),  some  of  the  classi- 
cal flow. examples  we  chose  have  a Re,\  increasing  with  x,  some  decreasing  with  x 
and  some  remaining  constant.  In  the  boundary  layer  Re,\  increases  in  the  down- 
stream direction  because  the  layer  thickness  6 ■ b increases  while  w remains 
constant.  The  same  is  true  for  the  free  shear  (mixing  layer)  and  for  the 
same  reasons.  The  same  is  also  true  for  the  two-dimensional  (2-D)  jet,  but 
only  because  b in  it  increases  faster  than  its  w decreases.  In  a second  class 
of  flows  of  constant  Re/y  belong  the  2-D  wake  and  the  axi- symmetric  (circular), 
or  A/S,  jet.  The  incompressible  A/S  wake  belongs  in  a class  by  itself  since 
its  Re,\  decreases  with  x,  i.e.,  from  the  body  on  downstream. 

If  compressibility  is  defined  as  the  case  where  the  temperature  decreases  in 
the  downstream  direction,  then  the  factor  (T/Te)"n  in  Equation  (5)  accelerates 
the  increase  of  Re/y  with  x,  or  causes  such  an  increase  in  the  flows  with  pre- 
viously constant  Re^.  For  the  2-D  wake,  for  example,  Re,y  will  now  increase  in 
the  downstream  direction,  since  the  wake  cooling  will  decrease  the  kinematic 
viscosity.  For  the  2-D  jet,  in  which  Re/y  increased  with  x even  when  T • Te 
because  of  the  rapid  spreading  of  the  jet,  the  cooling  of  the  gas  during  the 
expansion  will  accelerate  this  increase. 

The  effect  of  compressibility  on  Re/y  is  not  as  clear,  at  this  juncture,  for 
the  boundary  and  free -shear  layers;  although  their  increase  in  Re/y  with  x will 
not  be  affected  qualitatively,  the  rate  of  increase  will  be  affected  since  in 
these  cases  the  edge  conditions  do  not  change.  Separate  discussion  of  these 
cases  is  more  appropriate  when  the  details  of  these  flows  are  examined  later. 

» 

' The  axisynmetric  wake  remains  a special  - and  intriguing  - case  when  temperature 
effects  are  present.  Decaying  temperature  competes  with  the  decaying  wb  product; 
in  many  cases  an  increase  and  then  a decrease  of  Re/y  is  possible  here.  This 


case,  too,  will  be  discussed  later 
3.3  THE  TRANSITION  CRITERION 


Our  entire  approach  is  based  on  the  statement  that  no  turbulent  flow  is  possible 
if  Re/\  lies  below  a certain  magnitude  Re^o  which  is  universal.  Conversely,  we 
say  that  if  the  flow  is  turbulent  then  Re/y  > Re^0  everywhere  within  it.  The 
question  of  precisely  what  Re^0  is  will  be  discussed  below.  For  the  moment, 
let  us  see  what  this  statement  implies  for  the  comnon  turbulent  flow  examples 
described  just  above. 

Let  us  temporarily  return  to  an  example  with  an  obvious  history  of  Re^Q,  such 
as  that  of  the  incompressible  2-D  jet  idealized  so  that  it  can  be  represented 
by  Equation  (5)  with  T ■ Tg: 

b uebo  b 

Re.  ■ const  bw  Re'  - (const)  r—  w ( ) ■ C,  r—  w Re,  (6) 

A b0  Ve  4 bQ  bO 

In  this  expression  we  already  agreed  to  consider  C ^ a scaling  ratio  which  is  a 
pure  number  for  all  2-D  jets;  thus  the  Reynolds  number  Re^o,  where  b0  is  a 
constant  for  this  particular  jet  (e.g.,  the  jet  exit  size  or  the  jet  thrust 
integral)  controls  the  absolute  magnitude  of  Re,\.  The  variations 

b 


w ~ X/A 


(7) 

(8) 


— — /x 
bo 


(9) 


show  how  the  jet  turbulence  Reynolds  number  increases  from  the  jet  exit  on. 

Let  us  now,  with  the  aid  of  Figure  3,  visualize  the  Re;\  variation  of  such  a two- 
dimensional  turbulent  jet,  first  in  the  case  where  Rebo  is  so  large  that  Re/y 
lies  everywhere  above  Re,^0.  If  we  compute  Re^y  for  such  a given  jet  we  will 
find  Re/\  > Re^0  everywhere,  and  we  will  chus  find  no  reason  why  such  a jet 
cannot  exist.  Consider,  however,  the  jet  of  the  second  example  shown  in 
Figure  3,  where  Re/y  > Re/\0  beyond  a certain  distance  xQ  from  the  jet  origin, 
but  Re^  < Re.\0  for  x < x0.  If  such  a jet  had  been  postulated  as  everywhere 
turbulent,  this  postulate  would  have  been  invalid  for  a length  xQ  beyond  the 
jet  exit.  That  is,  if  we  began  checking  on  the  Re/y  of  this  jet  from  infinity 
forward  toward  the  origin,  we  would  arrive  at  a point  xQ  ahead  of  which  we 
would  find  the  molecular  viscosity  capable  of  dissipating  the  turbulence. 
Turbulence  is,  therefore,  impossible  ahead  of  xQ;  and  the  jet  would  have  to  be 
laminar  between  the  exit  and  x0. 


It  must  be  noted,  of  course,  that  the  condition  Ke\  > Rej\0  is  a necessary  one, 
but  not  a sufficient  one.  The  fact  that  the  jet  will  be  laminar  ahead  of  xQ 
is  established,  but  not  so  that  the  jet  will  be  turbulent  aft  of  x0.  For  that 
matter,  the  "high  Re^0"  example  of  Figure  3 may  also  be  laminar;  all  we  had 
said  previously  is  that  it  is  possible  for  this  case  to  have  a fully  turbulent 
jet  from  the  exit  on  to  infinity.  The  condition  of  sufficiency  is,  therefore, 
lacking  and  will  remain  so.  Nevertheless,  we  will  now  begin  referring  to  the 
physical  boundary  xQ  as  the  "transition  point,"  staying  aware  that  the  term 
refers  to  the  threshold  in  the  necessary  condition. 

The  arguments  built,  above,  around  the  2-D  jet  require  no  change  when  it  comes 
to  other  flows  of  the  same  class,  i.e.,  flows  of  monotonic  increase  of  Re^  with 
x.  These  flows  include  the  free  shear  layer  and  boundary  layer,  the  compressible 
A/S  jet  and  the  compressible  2-D  wake.  The  procedure  in  all  these  cases  is  the 
same:  we  first  assume  that  the  flow  is  fully  turbulent  and  then,  beginning  with 
the  far  end  of  the  flow,  we  calculate  the  Re^  proceeding  toward  the  flow  origin. 
We  may  never  find  chat  Re,  < Re/\o»  in  which  case  we  say  chat  it  is  possible  for 
the  entire  flow  to  be  turbulent  (or  we  may  risk  the  statement  that  the  entire 
flow  is_  turbulent).  Or,  if  we  come  Co  a point  where  Re,\  just  equals  ReAo,  we 
can  say  that  the  flow  ahead  of  that  point  is  laminar,  or  even  risk  the  statement 
Chat  transition  to  turbulence  occurs  at  that  point.  We  fully  intend  to  take 
the  latter  risk  in  this  work. 

The  situation  is  even  simpler  for  the  class  of  flows  with  constant  Re/V,  like  the 
incompressible  axi symmetric  jet  and  the  incompressible  2-D  wake.  In  these  cases, 
according  to  everything  said  so  far,  the  flow  is  either  wholly  laminar  or  wholly 
turbulent;  there  is  no  chance  that  transition  will  appear  at  some  distance 
downstream  of  the  flow  origin.  The  unique  case,  however,  of  the  axisymnecric 
wake  is  more  complex.  In  the  incompressible  version  of  this  flow,  Re,  decreases 
monotonically  and  two  possibilities  exist  according  to  our  Chinking:  "either 
chis  wake  is  wholly  laminar  or  it  is  turbulent  for  some  distance  behind  Che 
body  and  turns  laminar  beyond  some  "relaminarization , " or  reverse  transition, 
point.  In  the  compressible  case  (decaying  T in  Equation  (5)  ),  a maximum  in 
the  variation  Re,\(x)  may  even  exist  and  a turbulent  "piece"  of  wake  appears 
. somewhere  downstream,  preceded  and  followed  by  laminar  flow. 

Our  association  of  the  threshold  Re,  - Re,\0  with  transition  will  look  much  less 
premature  if  we  realize  chat  actual’  transition,  as  observed  to  occur  in  the 
sample  flows  discussed  above,  also  follows  the  same  rules.  For  example,  whereas 
2-D  incompressible  wakes  are  observed  to  be  either  wholly  laminar  or  wholly 
turbulent,  the  compressible  2-D  wake  can  become  turbulenc  at  some  downstream 
location.  Even  Che  odd  behavior  of  transition  in  a compressible  axisytanecric 
wake,  predicted  in  the  previous  paragraph  has  been  recognized  and  discussed  by 
Lees  (Reference  3).  Thus,  che  experimental  evidence  strongly  supports  our 
postulate,  even  beyond  its  stated  limitations.  What  we  have  here  is  the  ability 
to  make  a general  and  apparently  valid  statement  abouc  laminar- turbulent  tran- 
sition, which  is  based  on  physics  and  which  does  noc  require  "a  priori"  extensive 
• mathematical  development. 
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3.4  THE  NUMERICAL  VALUE  OF  THE  REYNOLDS  NUMBER  THRESHOLD 


To  make  quantitative  predictions  of  the  transition  occurrence,  we  need  a 
numerical  value  for  Re^0.  Such  a value  has  been  found  by  reversing  the  process, 
i.e.,  by  utilizing  transition  data,  in  this  instance  from  two-dimensional  wakes. 
The  derivation,  reported  in  Reference  9,  gave  a value  of  Re,\0  *■  15.  This 
magnitude,  consistent  with  the  intuitive  expectation  of  Re.\0  > > 1,  will  serve 
us  until  it  is  adjusted  (perhaps  only  slightly)  by  later  data  or  rationalized 
by  more  precise  mathematical  arguments. 

3.5  APPLICATIONS 

3.5.1  TRANSITION  IN  TWO  DIMENSIONAL  WAKES 

The  first  use  of  the  threshold  criterion  appeared  in  1978  in  Reference  9.  The 
application  involved  two  dimensional  wakes  for  which  a satisfactory  amount  of 
experimental  data  is  available.  The  results  are  outlined  below,  especially 
because  they  outline  the  general  method  of  application. 

The  point  of  departure  is  Equation  (1),  in  which  the  quantities  A,  v and  u'  are 
sought  to  be  expressed  in  terms  of  the  transition  discance  x-p  and  the  parameters 
of  the  wake  flow.  For  the  compressible  two  dimensional  wake  (Reference  9): 

u'(0)  ■ [ue  - u(0)]  (t  + 1)^  , C^  ■ constant  (10) 

where  (0)  refers  to  the  axis  (centerplane)  values  and  "e"  to  -he  edge  conditions, 
and  where  t is  the  "temperature  defect" 

(id 

*e 

For  the  scale  length  A we  obtain  from  Reference  9 


A 


C2  CDh 


ue  - u(0) 


(12) 


where  Cph  is  the  "wake  drag  thickness,"  i.e.,  the  integral  of  its  momentum  flux. 
The  kinematic  viscosity  in  (1)  is  furthermore  converted  to  temperature  via  the 
relation 


u 


(13) 


where  k is  fixed  by  the  thermodynamics  of  the  gas.  It  is  also  argued  that  at 
transition,  the  wake  (centerplane)  temperature  will  be  close  to  that  of  the 
laminar  wake  just  preceding  transition  and  if  the  latter  occurs  at  a distance 
x from  the  body,  then  laminar  wake  theory  gives 


/tj(V-l)(l+Q)Me2Rv 

4/Tt 


% 


t 


(14) 


where  a is  che  molecular  Prandtl  number,  V Che  specific  heac  ratio,  Mg  che  edge 
Mach  number,  Q che  non-dimensional  version  of  che  heac  Q exchanged  between  body 
and  flow 


Q * 


A 


- <3.ue3cDh 

2 


and  is  Che  wake  Reynolds  number 


R«w  * ^ CDh 


(15) 


(16) 


If  these  equations  are  inserted  into  (1)  ac  che  polnc  where  cransicion 
occurs,  i.e.,  when  Re.  " Re,0,  we  obtain 


Xp  a(v-l)2M.4(l-«3)2R«w 
C^h  " 16  tt 

'o 

R*o  " C& 


Re. 


1/k  + 0.5 
) -1 


-2 


(17) 


(18) 


Equation  (17)  is  che  desired  end  product;  it  is  a formula  giving  che  cransicion 

discance  xt  in  cerms  of  Che  given  parameters  of  Che  problem  k.  O’,  v,  M^,  Q, 

Rev  and  Cph;  Che  constants  and  C2  are  presumably  universal,  and  have  been 
determined  by  previous  experiments  with  turbulent  wakes  (e.g..  Reference  10). 
Equation  (17)  gives,  for  the  first  time,  che  dependence  of  che  transition 

distance  on  che  type  of  gas,  for  example,  X?  should  be  about  twice  as  large  in 

helium  (y  ■ 1.66)  chan  in  air  (y  ■ 1.4),  all  other  conditions  being  equal.  It 
also  explains  why  che  observed  transition  distances  increase  so  rapidly  with 
Mach  number  and  why  one  can  delay  transition  by  heating  che  wake- shedding 
body  (Q  > 0). 

Comparison  of  Equadon  (17)  with  experimental  data,  discussed  in  Reference  9, 
gave  good  results  and  convinced  this  author  that  further  work  on  che  two- 
dimensional  wake  is  not  needed. 

3.5.2  THE  TWO  DIMENSIONAL  BOUNDARY  LAYER 

Boundary  layers  were  addressed  after  che  work  of  Reference  9 was  completed, 
since  abundant  transition  data  for  comparison  exist  and  since  Che  boundary  layer 
cransicion  issue  is  of  great  current  interest.  This  study  is  continuing;  the 
following  text  describes  work  done  in  the  spring  and  summer  of  1978. 


To  apply  che  criterion  (1)  co  Che  boundary  layer,  calculation  of  che  quantities 
u',  A and  T(y)  (where  y is  the  usual  height-above-the-wall  coordinate)  is  again 
needed.  In  specifying  these,  existing  knowledge  is  utilized  ranging  from  known 
data  on  u'  and  A to  the  more  prosaic  textbook  formulas  for  u(y)  and  T(y)  in 
turbulent  boundary  layers.  Surprising  insight,  in  fact,  can  be  gained  by  mani- 
pulating commonplace  notions  on  those  boundary  layer  transition  features  fre- 
quantly  observed  in  experiments.  One  such  feature,  necessary  for  exploiting 
the  criterion  (1)  further,  is  che  geometry  of  che  transition  zone. 

3.5. 2.1  The  Virtual  Origin  of  the  Turbulent  Boundary  Layer  Flow.  We  can  make 
use  of  existing  knowledge  on  boundary  layers  to  compose  the  geometrical  picture 
of  transition  in  a boundary  layer  shown  in  Figure  4.  The  key  to  this  picture 
is  the  well-known  difference  in  boundary  layer  growth  rates  between  laminar 
and  turbulent  flows.  The  laminar  boundary  layer  edge  5 is  definable  despite 
che  diffuse  nature  of  the  laminar  interface  by  che  simple  expedient  of  picking 
a percentage  of  the  velocity  ratio,  e.g,  u(y*6)  ■ 0.99ue.  Similarly,  the 
turbulent  boundary  layer  edge  is  definable  despite  the  turbulent  intermittency , 
e.g.,  by  picking  6 to  correspond  to  a given  intermittency  factor. 

The  consequences  of  the  construction  of  Figure  4 are: 

(a)  A theoretical  transition  "point"  xj  is  defined  at  the  intersection  of 
the  laminar  and  turbulent  growth  curves. 

(b)  At  this  transition  "point"  the  laminar  and  turbulent  thicknesses  are 
che  same: 

5c(x  - xr)  ■ 6^  (x  - xT)  (19) 

(c)  The  virtual  origin  x of  the  turbulent  flow  lies  downstream  of  che  origin 
of  the  flow  (at  x ■ 0)  and  ahead  of  the  transition  point. 


Since  che  choice  of  ordinate  in  Figure  4 is  arbitrary,  we  can  replace  it  by 
the  integrals  6*  or  9: 


5*  (x  ■ xT) 
et  (x  - xj) 


5*  (x  » xx) 

(20) 

9-C  (x  - Xx) 

(21) 

so  that  the  shape  factors  such  as  (5/0)t  and  (6 /q)^  would  be  the  same  for 
laminar  and  turbulent  boundary  layers.  Calculations  done  in  the  present  con- 
tract period  for  typical  laminar  and  turbulent  boundary  layers  showed  that  this 
is  so  to  che  extent  chat  we  can  describe  che  turbulent  case  by  the  present 
empirical  approaches.  (See  Figure  5). 

To  act  as  further  reinforcement  of  che  construction  of  Figure  4,  we  have  cal- 
culated che  virtual  origin  of  che  turbulence  in  a boundary  layer  by  using  che 
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standard  form  of  the  momentum  thickness  (e.g..  Reference  11)  in  the  turbulent 
case, 

6 

, Re.  7 


9t  “ 12  (0-0262)  C (M) 


and  in  the  laminar  case, 

h ’ °-M 


where 


f<M>  3 < 7-1  2^ 


*T  " xc 


u.  ue  u 

Re ' ■ — , Re , ■ — A , Re„_  ■ — x-. 
v ’ A v * XT  v.  ” 


By  equating  Equation  (22)  with  (23),  according  to  Equation  (21),  we  obtain 


**  [RexX]5/12[f(M)]7/6 

The  results  plotted  in  Figure  6 show  that  at  low  Me  and  high  RexX,  the  virtual 
origin  x0  lies  very  close  to  the  transition  point  xx.  It  would  appear  that  in 
the  other  extreme  of  high  Me  and  low  RexX  the  virtual  origin  can  lie  upstream 
of  the  actual  flow  origin;  in  the  flat  plate  boundary  layer  (to  which  equations 
(22)  and  23)  in  anyway  refer)  this  would  mean  that  the  virtual  origin  lies  up- 
stream of  the  plat  leading  edge.  There  is  nothing  unreallistic  about  this 
concept,  since  the  virtual  origin  is  only  an  artifice  convenient  for  calculating 
the  flow,  and  there  are  many  situations  in  fluid  flows  where  the  virtual  origin 
lies  upstream  of  the  physical  generator  of  the  flow. 

Since  the  transition  Reynolds  number  Re^  is  known  to  depend  on  (in  a manner 
which  is  the  target  of  this  work),  the  actual  trace  of  (xx-x0)/xx  vs  Me  will 
cut  across  the  curves  of  Figure  6.  The  general  view  is  that  RexX  increases  with 
quite  rapidly  once  Me  > 5,  and  this  will  cause  this  trace  to  remain  below 
unity  in  the  figure.  An  example  of  this  is  shown  by  the  dashed  line  in  Figure  6 
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which  anticipates  the  forthcoming  results  of  RexT  vs  Me  for  the  flat  plate. 

It  is  seen  that  in  the  subsonic  and  supersonic  regimes  x0  lies  rather  close  to 
xj,  but  shifts  to  the  origin  of  the  flow  (i.e.,  xQ  becomes  nearly  zero)  at 
hypersonic  speeds. 

3. 5. 2. 2 Calculation  of  the  Reynolds  Number  Across  the  Flow.  We  now  return  to 
the  original  objective  of  applying  Equation  (1)  to  the  boundary  layer,  with  a 
view  to  deducing  the  transition  location  for  various  conditions. 

In  the  published  work  of  Reference  9,  the  Re/v  was  assumed  constant  across  a 
section  of  the  flow.  So  far,  this  idea  was  retained  in  the  previous  remarks 
because  we  wanted  to  put  emphasis  on  the  principle  and  to  discuss  various  shear 
flows  in  a general  way.  If  we  still  were  to  consider  Re/i  constant  across  a 
turbulent  boundary  layer,  we  could  quickly  derive  a transition  criterion  utilizing 
all  that  has  been  previously  said.  For  the  incompressible  2-D  layer  on  a flat 
plate,  for  example,  Equation  (10)  would  become  (since  u(0)  » 0): 

u'  - ClUe  (28) 

while  (12)  would  become,  using  "t"  to  denote  the  turbulent  case, 


A 


Vc 


(29) 


and  since  v ■ ve,  the  Re,\  at  any  position  along  the  layer  would  be 


ue5t 

R«j\  - cLc2  ( — ) 


(30) 


The  method  would  then  have  us  assume  chat  very  far  from  the  leading  edge  the 
layer  is  curbulenc  and  that  Re,  in  it  is  given  by  (1);  then  as  we  advance  toward 
the  leading  edge  Re,\  decreases  since  5C  decreases  until  we  arrive  at  Che 
threshold 


Re/\  - Re/v0  - 


(31) 


where  "T"  signifies  transition.  At  chat  point,  however,  5cT  - 5</p  according  to 
Equation  (19);  furthermore,  a more  familiar  form  is  obtained  by  using  the 
laminar  form  factor  C'  to  convert  to  6/T: 


i 


Re.Ao  - qC-C'  ( ) (32) 

e T 

The  quancicy  in  parenthesis ■ is  Che  familiar  momentum  Reynolds  number  ac  tran- 
sition based  on  the  laminar  momentum  thickness;  thus, 

*•.0 

Reof  ■ ■ constant  (33) 

C^C' 
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and  the  transition  momentum  Reynolds  number  for  the  incompressible  flat  plate 
boundary  layer  has  been  found.  Note  that  the  four  constant  quantities  forming 
Reqx,  although  empirical,  have  definite  physical  meaning  and  can  be  found  in 
the  literature  of  turbulence,  not  transition,  experiments. 

We  now  come  to  grips  with  the  fact  that  Re  is  not  constant  across  the  flow, 
but  varies  with  y as  u',  A and  v vary  withvy*  Suitable  expressions  for  u'(y), 
etc.,  can  be  found  in  the  present  bibliography  on  turbulent  boundary  layers, 
which  is  considerable.  The  following  describes  the  work  done  so  far  to  pick 
these  functions. 

(a)  Velocity  Fluctuations 

The  u'(y)  variation  was  obtained  by  examining  experimental  data  available  in 
the  literature.  Guidance  for  the  most  rational  correlation  of  these  data  was 
provided  by  the  original  work  of  Morkovin  (Reference  12)  according  to  which 
u'  should  scale  with  the  friction  velocity  and  a compressibility  factor.  The 
rationale,  together  with  a list  of  data  sources  is  given  in  Reference  13.  In 
that  reference  it  is,  therefore,  pointed  out  that  the  correlated  form  is 


ci  > 


TW 

p(y) 


(34) 


Inspection  of  the  data  revealed  that  the  function  is  additionally  sensitive 
to  the  wall-to-stagnation  temperature  ratio  TW/TQ  but  independent  of  the 
Reynolds  number.  The  function  was  curve- fitted  as  follows: 

T 2 * 3 

C1  < 7 . f2  > " Cl. 814-0. 2188516  j - 2.797874  ( r ) + 1.430676(  |)  ] 


x (0.4  + 0.6  — ) 
1o 


(35) 


(b)  The  Longitudinal  Scale  of  Turbulence 

Information  of  scale  lengths  is  extremely  meager;  on  the  basis  of  this  author's 
earlier  experiences,  the  following  was  adopted,  Independent  of  y/5: 

A - 0.26 t - C25c  (36) 

(c)  Kinematic  Viscosity 

pe  , T " 

v “ '■’e  p 4e  ; (37) 

where  we  have  assumed,  quite  safely,  that  the  layer  is  isobaric.  The  conscant 
n is  taken  to  be  0.75  in  subsequent  computations. 
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(d)  Other  Variables 

When  re-casting  Equation  (32)  into  a form  suitable  for  computations,  we  need 
the  following  additional  properties  of  the  turbulent  boundary  layer:  velocity 
profile  form,  temperature-velocity  relations,  skin  friction  variation,  and  a 
connection  between  skin  friction  and  momentum  thickness.  Unlike  the  turbulence 
properties  u'  and  A,  these  additional  properties  are  supported  by  a large 
volume  of  experimental  data.  Improvements  and  additions  of  such  data  are  con- 
stantly appearing  in  the  literature,  and  the  present  analysis  could  be  cer- 
tainly postponed  until  these  properties  are  determined  with  great  precision  for 
a wide  variety  of  boundary  layers.  In  order  to  proceed,  therefore,  we  limited 
our  approach  to  zero-pressure-gradient  flat  plate  boundary  layers  in  0 < Me  < 10, 
0.2  < tw/TQ  < 1 in  exchange  for  using  established  formulas  for  the  needed  para- 
meters as  follows: 


- for  the  skin  friction, 


2 Tv  1 

Cf  - j - 0.0262  f (Me) 

P u 
e e 


6/  ReT77 


with 


f(Mj  3 ( 


5/7 


e'  Y-l 


) 


(38) 


(39) 


uex 

Re„  3 

x ve 


(40) 


- for  the  momentum  thickness  growth, 


d9 


dx  P u 

e e 


(41) 


which,  with  the  aid  of  (38),  integrates  to 


(42) 


- for  the  velocity  profile, 
1/7 

~ ■ (y/5) 


(43) 
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- for  the  temperature-velocity  relation,  the  Crocco  relation 


T Tw  Y-l  j u Y”  1 o u 2 

+(1+  — Me-^~)  — - — Me(  — ) (44) 

Using  Equations  (38),  (41),  (43)  and  the  isobaric  assumption  p = pe  into 
Equation  (32),  we  now  obtain 


Vt  cf  p n+1 

ReA  " C2C1  TT  ) v T~  ( T"  ) 


M 2 ' pe 


Since  it  is  customary  to  refer  Reynolds  numbers  to  the  momentum  thickness, 
the  following  conversions  can  be  next  made  using  the  formulas  previously 
presented: 

7 7 l -l 

Cf  - C36  Cf(Me)]6  ( -jj)  * Re9  6 5 c3  " 0-0262  (46) 


ue9t 
Re.  = 

9 ve 


Ve  " 9t  R69 


so  that  the  final  product  is: 


7-11/12  " °*°1767  C1  “ < “ > 


Thus,  we  see  that  for  a given  turbulent  boundary  layer,  the  transverse 
variation  of  Re/^  depends  on  the  corresponding  transverse  variations  of  the 
fluctuations  via  (given  in  Equation  (35)  ) and  of  the  density  ratio  ob- 
tained by  combining  (43)  with  (44).  Its  edge  Mach  number  dependence  enters 
explicitly  via  f(Me)  and  implicitly  via  p/Pe  and  the  shape  factor  5/9.  The 
Re,  depends  on  the  ratio  TW/TQ  through  its  effect  on  Cj^  (Equation  (35)  ) and 
°/p»  (by  the  Crocco  relation).  Finally,  the  dependence  on  Reg  is  almost 
wholly  left  to  the  factor  dividing  ReA  in  Equation  (49);  the  dependence  of 
5/9  on  Rea  is  apparently  negligible. 
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3. 5. 2. 3 Computations  of  the  Turbulence  Reynolds  Number.  Using  Equation  (49), 

the  ratio  on  the  l.h.s.  was  computed  for  Mfi  * 0,  1,  2 9,  10  and  TW/TQ  = 

0.2,  0.4,  0.6,  0.8  and  1.  Typical  results  are  shown  in  Figures  7 through  9. 

On  the  wall,  the  value  of  Re«,  (or  Re^/Reg^^ is  controlled  primarily  by  the 
wall  temperature,  while  its  value  at  the  layer  edge  is  controlled  by  the  mag- 
nitude of  the  fluctuations  there.  Since  the  fluctuations  (i.e.,  C^)  decrease, 
and  the  density  ratio  increases,  toward  the  edge  Re\  usually  has  a maximum 
value  at  some  point  intermediate  to  y =*  0 and  to  y * 5 . Also,  while  Re/\  is 
nearly  proportional  to  Reg,  it  decreases  as  Mg  is  increased  and  also  as  the 
wall  is  cooled.  For  example,  in  the  adiabatic  case 


Maximum 




11/12 


= 0.33  at  M_  =* 


=*  0.09  at  Mg  =*  10  (51) 

Similarly,  at  M„  *■  10  the  maximum  Re,\  value  decreases  from  0.09  in  the  adiabatic 
case  to  about  0.041  at  Tw/Tq  = 0.2. 

3. 5. 2. 4 Transition  Predictions  from  the  Computed  Reynolds  Numbers.  In  contrast 
to  the  hypothetical  example  of  Section  3. 5. 2. 2 where  Re,\  was  assumed  independent 
of  y,  the  calculations  described  above  show  that  Re.\  can  vary  considerably  in 
the  range  0 < y < o.  The  meaning  of  this  variation  is  necessarily  complex. 

The  rules  established  can  be  explained  with  the  aid  of  Figure  10.  The  boundary 
layer  in  question  is  first  assumed  turbulent  throughout  so  that  its  variation 
Reg(x)  is  known.  Thus,  at  any  x,  a curve  Re/\(y/5)  can  be  plotted,  like  those 
of  Figures  7 through  9,  except  that  the  ordinate  is  replaced  by  just  Re,\.  As 
done  before,  this  Re,\  is  compared  with  the  threshold  Ret0,  and  three  possibilities 
exist : 

(a)  Re;\0  > Rei\  everywhere  (case  (a)  in  Figure  10).  The  boundary  layer  is 
laminar  at  the  position  x where  this  Re,\  was  computed. 

(b)  ReA0  < ReA  everywhere  (case  (b)  in  Figure  10).  The  boundary  layer  can  be 
turbulent  at  this  x;  we  assume  that  it  is. 

(c)  The  Re,\0  and  ReA  curves  cross  (case  (c)  in  Figure  10).  Strictly  speaking, 
no  statement  can  be  made  in  this  case  based  on  the  present  approach,  since  the 
Re,\  curve  in  this  case  was  computed  (from  Equation  (49)  ) assuming  that  the 
layer  was  turbulent  for  all  y/5.  However,  we  can  see  that  if  the  layer  was 
turbulent  everywhere  in  the  sense  that  Equations  (38),  etc.,  are  valid,  then 
Re,\  would  indeed  have  a physically  possible  value,  which  is  larger  than  Re^g. 

In  other  words,  there  is  no  turbulent  boundary  layer  which  could  produce  such 
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a case  as  shown  by  (c)  in  Figure  10.  This  case,  Coo,  must  then  represent  a 
"non- turbulent"  boundary  layer. 

The  above  considerations  led  us  to  the  only  remaining  possibility,  i.e.,  that 
transition  (departure  from  turbulence)  occurs  as  soon  as  the  Re \ line  first 
touches  the  Re/\  curve  (case  (d)  in  Figure  10).  But  if  this  is  our  definition 
of  transition  chen,  Co  be  consistent,  the  Re0  value  (of  this  specific  boundary 
layer)  where  this  happens  must  be  that  based  on  the  laminar  9,  according  to 
Figure  4.  Thus,  the  value  of  ReQ  at  which  this  first  coincidence  of  the  two 
curves  occurs,  is  the  standard  laminar  Reg  of  transition. 

Following  this  definition,  the  transition  Reynolds  numbers  Rea-p  were  computed 
for  boundary  layers  in  the  chosen  range  of  0 < Mg  < 10,  0.2  <TW/T0  < 1.  The 
results,  shown  in  Figures  11  through  13,  display  some  quantitative  and  some 
qualitative  features  which  are  remarkable  for  their  similarity  to  common 
observations.  These  features  are: 

1.  For  fixed  Tw/Tg,  Regx  first  decreases  with  Me  and  then  increases.  The 
edge  layer  (y/6  = 1)  is  responsible  for  the  decreasing  portion  and  che  wall 
layer  (y  “•  0)  for  Che  increasing  portion. 

2.  Below  Me  *•  4.5  cooling  always  delays  transition  (increases  RePx,  "stabilizes 
the  layer").  Above  Mg  *■  4.5  transition  reversal  occurs.  The  reversal  Tw/Tq 
decreases  as  Me  increases;  at  large  Me  (of  order  10  for  example),  cooling 
moves  transition  forward  (to  smaller  Regx)  until  a rather  low  Tw/T0  is  reached, 
below  which  Tw/T0  has  a stabilizing  influence  again. 

3.  It  follows  from  the  above  that  the  reversal  phenomenon  (9RePx/bMe  or 
9RegT/S(Tw/T0) changes  sign)  occurs  whenever  the  transition-ending  role  is 
switched  from  one  y/5  position  in  the  layer  to  another  (e.g.,  from  the  wall 
layer  to  che  edge  layer  or  vice  versa). 

4.  The  transition  scaling  parameters  are  largely  the  same  as  those  scaling  the 
turbulent  boundary  layer,  i.e.,  those  factors  affecting  che  friction  coefficient 
Cf,  profile  b/pg  and  distribution  C^(y/S)  of  the  fluctuations  in  a turbulent 
layer,  etc.,  also  affect  transition. 

3.5.3  DISCUSSION  OF  THE  BOUNDARY  LAYER  TRANSITION  RESULTS 

The  results  shown  in  Figures  11  through  13  were  obtained  by  postulating  chat 
transition  occurs  when  che  straight  Re y0  lines  on  Figure  10  first  encounter  che 
bottom  part  of  the  Re/\  curve  obtained  at  chat  position  x in  the  boundary  layer. 
Any  further  increase  in  Reg  will  raise  the  Re(\  curve  above  the  Re.\0  curve,  and 
the  situation  will  chen  resemble  case  (a).  Since  the  latter  case  represents 
the  "established"  turbulent  boundary  layer,  it  is  clear  Chat  if  our  transition 
criterion  is  valid,  it  is  a criterion  for  the  "end"  or  completion  of  transition 
rather  chan  for  its  "beginning"  or  "onset".  So  we  say  that  every  part  of  che 
layer,  from  y ■ 0 to  6 , will  have  to  meet  the  condition  Re,\  £ Re.\0  before  the 
entire  boundary  layer  is  considered  turbulent. 
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Ic  is  Chen  cempcing  Co  consider  che  case  (e)  of  Figure  10,  where  Rep  is  jusc 
large  enough  Co  bring  Che  Re,\  curve  inCo  firsc  conCacC  wich  Che  Re/\0  limic,  a* 
che  "beginning"  of  cransicion.  This  suggescion  cannoC  be  caken  coo  seriously 
because  of  Che  argumenCs  in  SecCion  3. 5. 2. 4;  however,  one  credible  feacure  of 
Chis  suggescion  is  ChaC  che  zone  of  che  boundary  layer  aC  y = y (Figure  10) 
is  che  place  in  che  layer  where  Che  firsc  appearance  of  Che  turbulence  will 
occur.  Thac  Curbulence  will,  in  face,  firsc  appear  in  specific  zones  of  Che 
boundary  layer,  racher  Chan  uniformly  across  ic,  a phenomenon  which  has  in- 
trigued observers  for  years  (Reference  14) . Figure  14  plots  the  location  of 
y0/5  from  Che  present  calculations.  This  picture  shows  not  only  chac  y 0/6 
increases  wich  Mg  and  Tw/T0  has  has  been  observed,  buC  also  Chat  a certain 
"switching",  triggered  by  heat  cransfer,  occurs  at  low  Mg.  The  increasing 
y0/5  as  Me  increases  is  also  in  general  agreement  wich  the  "critical  layer" 
idea  of  hydrodynamic  scabilicy  theory,  which  is  also  plotted  in  Figure  14. 

In  addition  to  the  work  necessary  co  arrive  at  Figures  11  through  13,  several 
issues  related  to  it  were  examined.  One  such  issue  refers  Co  che  shape  of  the 
velocicy  profile  and  ics  effect  on  the  density  profile.  Clearly,  Equation  (43) 
i.e.,  che  sevench-power  law,  is  che  most  elementary  one  which  can  be  assumed, 
and  ic  can  be  improved  by  replacing  it  with  the  Clauser-Coles  logarichmic 
velocity  formulation.  Also,  the  Crocco  relation  is  equally  crude  and  can  be 
replaced  for  example,  by  che  Whitfield-High  relation  (see  Reference  15).  Such 
changes  were  attempted  without  any  major  changes  in  the  results  of  Figures  11 
through  13,  so  far. 

A matter  of  concern  is  the  precise  values  of  che  quantities  on  the  r.h.s.  of 
Equation  (49).  For  example,  the  variation  of  u'  across  che  layer,  as  given  by 
the  curve- fit  of  Equation  (35),  ignores  che  existence  of  che  laminar  sublayer. 
At  the  boundary  layer  edge,  also,  the  value  of  che  fluctuations  u'  is  of  order 

0. 2.  since  the  edge  values  are  critical  in  deducing  Rea£,  improving  the 

u'(y  = 5)  should  have  some  effect  on  che  ReoT  deduced  so  far. 

3.6  PRESENT  ISSUES,  OBJECTIVES  AND  ACTIVITIES 

There  is  little  doubt  thac  the  simple  physical  criterion  as  expressed  by 
Equation  (1)  has  already  been  very  successful  in  making  valid  and  practically 
useful  predictions  of  transition  onset.  This  is  not  only  attested  to  by  the 
wake  resulcs  (Reference  9)  but  also  by  the  boundary  layer  results  of  Figures 
11,  12  and  13.  The  present  and  future  work  to  be  done  on  this  contract  de- 
serves considerable  thought  because  of  this  very  success. 

In  view  of  the  necessarily  restricted  contract  resources  and  the  requirements 
for  application  arising  from  various  engineering  centers  of  the  Air  Force,  the 
following  three  areas  are  considered  important  co  pursue: 

1.  Strengthening  of  che  underlying  concepts  of  the  theory. 

2.  Extension  to  boundary  layers  wich  other  transition-causing  factors  beside 
Me  and  Tw/T0,  e.g.,  freestream  Curbulence,  pressure  gradient,  etc. 
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3.  Application  to  flows  other  than  boundary  layers,  e.g.,  the  free  sheax 
layer. 

We  are  presently  pursuing  the  first  item  of  this  list.  The  sublayer  is  being 
taken  into  account  by  modifying  the  curve-fit  of  Equation  (35)  to  account  for 
the  sublayer  existence  and  thickness.  The  fact  that  the  fluctuations  do  not 
go  to  zero  at  y = 5 is  also  of  concern;  the  fact  itself  is  not  in  doubt,  but 
it  must  be  rationalized  on  valid  physical  grounds.  The  problem  here  is  analo- 
gous to  the  differences  in  lateral  spreading  of  heat  and  momentum  in  a boundary 
layer,  i.e.,  a Prandtl  number  must  exist  which  controls  the  lateral  spread  of 
turbulence  as  well  as  of  heat. 

In  the  near  future,  we  will  attempt  to  extend  Equation  (49)  to  account  for 
pressure  gradient,  surface  roughness  and  freestream  turbulence.  It  is  clear 
that  these  will  affect  ReQj  by  affecting  separately  the  skin  friction  and 
fluctuation  intensity.  Preliminary  work  has  shown  that  the  present  con- 
cept is  perfectly  capable  of  explaining  the  acceleration  of  transition  by  these 
three  agents.  The  key  to  success  here  is  the  existence  or  non-existence  of 
turbulent  boundary  layer  data.  That  is,  to  find  how  surface  roughness  affects 
transition,  it  is  necessary  to  know  the  roughness  effect  on  Cf  for  turbulent 
boundary  layers. 

Also  on  the  schedule  is  an  attempt  to  apply  the  present  approach  to  the  free 
shear  layer,  since  this  problem  is  of  interest  to  designers  of  systems  contain- 
ing heterogeneous  mixing.  The  chief  obstacle  here  is  the  lack  of  shear  layer 
turbulence  data  for  use  in  the  theory,  as  well  as  of  shear  layer  transition 
data  for  comparison  with  the  forthcoming  results. 

4.  HYDRODYNAMIC  STABILITY  OF  THE  HYPERSONIC  BOUNDARY  LAYER 

In  the  present  period  the  final  analysis  was  done  and  the  results  published, 
of  the  hypersonic  boundary  layer  experiment  performed  in  1976-77.  This  work 
represents  a complementary,  and  more  conventional,  approach  to  the  boundary 
layer  transition  problem  than  adopted  in  the  "dissipation  theory"  project 
described  in  Section  3 . Whereas  the  latter  is  a theoretical  approach  to 
transition  dealing  with  a necessary  condition  for  the  existence  of  turbulence, 
che  stability  experiment  gathered  data  dealing  with  the  mechanism  of  transition 
onset . 

Stability  experiments  under  AFOSR  sponsorship  at  Aeronutronic  began  in  1975, 
and  have  already  produced  nearly  half  the  existing  information  on  hypersonic 
boundary  layer  disturbance  growth  and  amplification.  The  program  takes  ad- 
vantage of  che  author’s  graduate  research  in  the  stability  of  laminar  boundary 
layers  on  flat  plates  at  Me  » 5.8  (Reference  16).  In  Che  1975-76  period,  this 
investigator  analyzed  and  published  (Reference  17)  data  taken  in  Che  laminar 
and  transitional  flow  over  a 5°  half-angle  cone  taken  at  M,  > 8 in  che  AEDC 
Wind  Tunnel  B.  The  two  experiments,  referenced  just  above,  along  with  the 
work  of  Kendall  (Reference  18)  hinted  at  a very  complex  stability  diagram  at 
hypersonic  speeds. 
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The  present  work  was  undertaken  partly  in  order  to  expand  the  Reynolds  number 

and  frequency  ranges  covered  in  References  17  and  18,  and  partly  in  order  to 

discern  the  effect  of  wall  cooling  on  stability.  The  experiment  was  done  in 
the  spring  of  1977  in  Wind  Tunnel  B at  AEDC,  using  a 4°  half-angle  cone  at  two 
temperature  ratios,  Tw/T0  ”0.4  and  0.8.  The  results  were  analyzed  in  the 
present  contract  period  and  presented  at  the  1978  Heat  Transfer  and  Fluid 
Mechanics  Institute  in  June  1978  (Reference  19).  Briefly,  the  results  were 
as  follows. 

1.  The  data  obtained  in  this  experiment  are  consistent  with  all  other  data 

reported  to  date  at  ^ 7.  A complete  picture  of  disturbance  amplification 
is  thus  available  at  this  M®  in  the  range  1100  < R < 2400  (or  1.2  x 106  < Rex 

< 5.8  x 10°)  and  for  TW/TQ  down  to  0.4.  Moreover,  the  consistency  of  the  re- 

sults of  three  independent  experiments  done  with  essentially  different  wind 
tunnels  and  models,  removes  suspicion  that  the  instability  was  grossly  in- 
fluenced by  peculiarities  of  the  flow  or  the  test  article. 

2.  The  present  results  fully  vindicate  the  admonition  of  Mack  and  others 
that  boundary  layer  stability  at  hypersonic  speeds  is  a complex  phenomenon. 

No  fewer  chan  three  distinct  unstable  regions  have  been  found  in  contrast  to 
the  single  instability  "loop"  encountered  at  subsonic  or  supersonic  speeds. 
Unstable  disturbances  exist  for  F > 4 x lO-^,  or  considerably  higher  in  fre- 
quency chan  previously  thought. 

3.  The  boundaries  and  amplification  race  isotherms  have  been  isolated  with 
considerable  precision.  The  "lower  unstable  region"  has  been  mapped  over  the 
entire  range  of  1100  < R < 2400.  A lower  neutral  branch  for  this  region  has 
been  tentatively  identified  but  its  minimum  Reynolds  number  has  not  been 
found.  By  inference  from  the  work  of  Mack,  this  region  is  the  one  generated 
by  inclined  first  mode,  and  normal  second  mode  instabilities.  It  is  clear 
chat  the  first  mode,  which  dominates  instabilities  at  lower  Mach  numbers,  is 
rather  unimportant. 

4.  The  middle  unstable  region,  newly  mapped  in  this  work,  serves  to  amplify 
disturbances  at  considerably  higher  frequencies  than  the  first  region.  What 
is  surprising  about  this  region  is  its  rather  large  amplification  rates; 
these  are  nearly  Che  same  as  chose  found  in  the  lower  unstable  region,  i.e., 
of  order  2 to  5 x 10“^.  Precisely  what  instability  mode  is  active  here  is 
presently  unknown. 

5.  A third  unstable  region  (here  called  "upper  unstable")  also  exists. 
Information  of  this  region  is  still  fragmentary  but  it  must  involve  wavelengths 
of  order  0.26  or  smaller. 

6.  The  passage  of  each  Fourier  component  through  Che  unstable  regions  results 
in  the  formation  of  Che  well-known  "laminar  waves."  Of  these,  the  prominent 
ones  belong  to  the  lower  unstable  region  (and  to  the  second  instability  mode) 
while  Che  growth  of  a weaker  "second  harmonic"  has  also  been  observed,  assoc- 
iated with  the  middle  unstable  region.  Physically,  the  waves  are  about  twice 
as  long  as  the  boundary  layer  is  thick.  The  significant  feature  of  these 
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waves  is  Che  new  finding  Chat  the  locus  of  their  maximum  amplitude  on  the 
stability  diagram  eventually  crosses  the  upper  neutral  branch  of  the  lower  in- 
stability region.  By  the  definition  of  the  neutral  branch,  this  crossing 
signals  the  end  of  the  growing  period  of  the  laminar  waves  and  the  onset  of 
boundary  layer  transition. 

7.  The  act  of  lowering  the  surface  temperature  from  0.8  to  0.4  has  only  a 
negligible  effect  on  the  location  and  size  of  the  unstable  regions  and  a 
similarly  slight  effect  on  the  laminar  waves,  but  the  amplification  rates 
nearly  double  in  both  the  lower  and  the  middle  unstable  regions.  This  overall 
increased  rate  of  amplification  could  only  mean  earlier  transition  to  turbu- 
lence, exactly  as  observed.  This  marks  the  first  occasion  where  the  cooling 
effects  on  amplification  and  transition  were  recorded  simultaneously;  their 
mutually  consistent  behavior  should  renew  confidence  in  the  relevance  of 
hydrodyuamic  stability  to  the  boundary  layer  transition  problem. 

5.  PUBLICATIONS  AND  PRESENTATIONS 

In  Che  period  October  1977  through  September  1978,  the  following  publications 
were  prepared  on  the  work  performed  under  this  contract. 

5 . 1 TECHNICAL  REPORTS 

- Demetriades,  A.,  "Experimental  Test  of  the  Theory  of  Multiple-Nozzle 

Cusp  Wakes,"  Aeronutronic  Publication  U-6395,  Ford  Aerospace  & Communications 
Corp.,  Newport  Beach,  California,  December  1977  (also  to  appear  as  AFWL  TR) . 

- Demecriades,  A.,  "Laminar  Boundary  LayeT  Stability  Measurements  at  Mach  7 
Including  Wall  Temperature  Effects,"  AFOSR  TR  No.  77-1311  (Aeronutronic 
Publication  No.  U-6381),  November  1977. 

5.2  JOURNAL  ARTICLES 

- Demecriades,  A.,  "New  Experiments  on  Hypersonic  Boundary  Layer  Stabisity 
Including  Wall  Temperature  Effects,"  Proceedings  of  the  1978  Heat  Transfer 
and  Fluid  Mechanics  Institute  (Pullman,  Washington,  June  26-28,  1978), 
Stanford  University  Press,  1978,  p.  39. 

- Demetriades,  A.,  'Transition  to  Turbulence  in  Two-Dimensional  Wakes," 

AIAA  J.  Vol.  16,  No.  6,  June  1978,  p.  587. 

- Demetriades,  A.,  "Turbulent  Fluctuations  in  Che  Wakes  of  Gas- Dynamic -Laser 
Nozzle  Cusps,"  J.  of  Energy,  Vol.  2,  1978. 

5.3  PRESENTATIONS  AT  MEETINGS  AND  SYMPOSIA 

- Demecriades,  A.,  'ftew  Experiments  on  Hypersonic  Boundary  Layer  Stabisity 
Including  Wall  Temperature  Effects,"  presented  at  the  1978  Heat  Transfer  and 
Fluid  Mechanics  Institute,  University  of  Washington,  Pullman,  June  26,  1978. 
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• Demetriades,  A.,  "Basic  Research  in  High  Speed  Aerodynamics,"  delivered  at 
seminar.  Mechanical  Engineering  Department,  Montana  State  University,  at 
Bozeman,  October  27,  1978. 

- Demetriades,  A.,  "Consequences  of  Imposing  a Necessary  Condition  for  the 
Appearance  of  Transition,"  presented  at  the  31st  Annual  Meeting,  Division 
of  Fluid  Dynamics,  University  of  Southern  California,  Nov  20,  1978. 

Interaction  with  Government  Laboratories 


- Space  and  Missile  Systems  Organization  (SAMSO),  LAAFS:  Exchange  of  personnel, 
consultations  and  technical  direction  meetings,  chiefly  on  related  Contract 
F04701077-C-0113.  The  dissipation  theory  of  transition  has  been  presented 

as  a candidate  scheme  for  explaining  blunt  body  boundary  layer  transition 
data  obtained  in  the  above  contract. 

- Arnold  Engineering  Development  Center  (AEDC),  Tullahoma,  Tennessee.  Exchange 
of  visits  relating  to  boundary  layer  transition  (J.  Donaldson,  S.  Pate, 

S.  Dougherty). 

- Air  Force  Weapons  Laboratory  (AFWL) , Albuquerque,  N.M.  (P.J.  Ortwerth). 
Exchange  of  visits,  consultations  and  technical  direction  regarding  laser 
nozzle  cusp  wake  research  and  research  on  wake  and  shear  layer  transition. 

- Wright  Air  Development  Center,  Aeronautical  Research  Laboratory  (WADC/ARL) 

(W.  Hankey,  J.  Shang).  Fluctuation  and  other  data  obtained  in  pressure 
gradients  in  Contract  F33615-77-C.-3016  are  useful  as  inputs  to  the  dissipation 
theory  of  transition. 
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Figure  2.  Typical  experimental  results  from  the  GDL  multiple  nozzle 

wake  study,  in  this  case  of  the  velocity  defect  (from  Ref.  6). 
The  theory  was  also  developed  in  the  present  research  (from 
Ref.  1). 
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Figure  3 


Qualitative  variation  of  turbulence  Reynolds  number  along  a 
two-dimensional  turbulent  Jet,  showing  how  the  turbulence- 
onset  criterion  applies  in  the  transition-prediction  theory 
developed. 


Figure  5.  Form  factor  computed  for  adiabatic  turbulent  boundary  layers. 
The  integrated  results  are  closer  to  the  experimental 
evidence  and  are  also  identical  with  theoretical  calculations 
for  laminar  flow. 


Figure  6.  Location  of  virtual  origin  of  boundary  layer  turbulence 

relative  to  the  transition  point  location,  for  flat-plate 
adiabatic  boundary  layers. 
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Figure  12.  The  results  of  Figure  11  as  recast  in  terms  of  wetted  length 
transition  Reynolds  number  instead  of  Reg. 
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